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ABSTRACT 

Accurate theoretical colllslonal excitation rates are used to determine 
the emissivities of CO rotational lines for 10 4 cm -3 < n(H2), 100 K < T < 2000 K, 
and J < 50. An approximate analytic expression for the emissivities which is 
valid over most of this region is obtained. Population inversions in the lower 
rotational levels occur for densities n(H 2 ) ** 10 3 ” 5 * cm” 3 and temperatures 
T > 50 K. Interstellar shocks observed edge-on are a potential source of 
millimeter-wave CO maser emission. The CC rotational cooling function suggested 
by Hollenbach and McKee (1979) is verified, and accurate numerical values given. 
Application of these results to other linear molecules should be straightforward. 
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I. INTRODUCTION 

The development of far-infrared spectroscopy is opening a new window on 
the Universe. Rotational lines of astronomically abundant molecules often lie 
in the far-infrared. Observation of these lines provides a powerful probe for 
determining both the density and temperature of the emitting gas, as recently 
demonstrated by the interpretation of the observations of the J ■ 21 ♦ 20, 

22 ♦ ti, 27 ♦ 26, and 30 ♦ 29 lines of hot «.0 in Orion (Watson et al 1980, 
Storey et al 1981). Using collislonal excitation cross-sections calculated by 
one of us (S.G.), Storey et al (1981) calculated the CO rotational level 
populations and line intensities for a number of temperatures and densities 
and showed that the observed intensities in Orion could be modelled by 
emission from two components: a 2000 K component with n(H 2 ) H 1 x 10^ cm"^ 
and a 400 - 1000 K component with n(H 2 ) 2 5-2 x 10^ cm"3. 

Inverting a 51 by 51 matrix to determine the level populations (as done by 
Storey et al) is somewhat cumbersome for interpreting data, and is completely 
Impractical for theoretical calculations of the emission spectra of shocks, 
etc. In this paper accurate numerical results for the line intensities for 
optically thin CO collisionally excited by H 2 are presented for a number of 
points in the regime n(H 2 ) > 10 4 cm -3 and 100 K < T < 2000 K (§11). In §111 we 
use an approximation Introduced by Hollenbach and McKee (1979; hereafter, HM) 
to obtain an analytic expression for the level populations and line intensities 
which is accurate to within a factor 1.3 over most of this regime. Since the 
theoretical collision cross-sections have an accuracy of about 30Z and 
astronomical observations average over regions at a variety of temperatures and 
densities, the analytic approximation should be quite adequate for astronomical 
purposes. Inversion of the lower rotational levels is discussed in §IV. We 
also use the numerical results to evaluate the 
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accuracy of the CO cooling rate found by HM (S V). 

Although we have focuaaed on CO as the typically most abundant molecule 
other than H 2 , our methods should also apply directly to other linear mole- 
cules such as HCN, HC 3 N, and HC0 + . 

II. EXACT CALCULATION OF POPULATIONS 

Calculation of the equilibrium population of each state requires first of 
all a detailed knowledge of the entire matrix of collislonal rate coefficients, 
Yjj', from every state, J, to every other state, J'. Although such detailed 
Information cannot yet be obtained experimentally, It has been possible to 
obtain fairly accurate values theoretically. This Involves two steps, computing 
the lntermolecular forces and then calculating the scattering dynamics. In an 
early study, Green and Thaddeus (1976) obtained the forces from a simple quantum 
mechanical model and solved for the dynamics by an essentially exact numerical 
method which Is tractable only at very low collision energies where few rota- 
tional levels are accessible. Subsequently Green and Chapman (1978) extended 
these results somewhat by using approximate scattering techniques. In a more 
recent study, Thomas et al. (1980) have recomputed the lntermolecular forcer via 
accurate ab lnltlo solution of the Schroedlnger equation; these authors showed 
that the improved potential predicted atate-to-state rates which differed In 
some cases by up to a factor of two from earlier values, although the total 
excitation rate did not change very much. As discussed by Green and Thomas 
(1980) available pressure broadening data can be used to provide some check on 
the theory, and. In fact, agreement Is within experimental error. 

At high energies, where many levels are accessible, the problem becomes 
unwieldy due to the number of state-to-state rates which must be considered. 
Fortunately, when the kinetic energy Is large compared to the rotational 
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energy spaclngs, It can be shown chat Che various races are not all 

Independent; In particular the entire matrix of rates can be obtained simply 

In terms of one row yjo (Goldflam et al. 1977): 

Y tt , - (2J'+1) l <2L+1)( J L J ') 2 Y , J’ < J (2.1) 

JJ L 000 LO 

where (J L J' ) is the usual three- j symbol. Eq. (2.1) applies only to 
^0 0 0 ' 

downward rates, but upward rates can be obtained from detailed balance, 

Yj f j - Y JJ| [(2J+1)/(2J , +1)] exp(-AE/kT) (2.2) 

For a given collision energy, this approximation becomes less accurate for 
higher rotational levels since the energy splittings increase with Increasing 
J. As shown by DePrlsto et al. (1970) one can approximately correct for this 
defect by a simple modification of Eq. (2.1) which takes Into account the 
rotational energy splittings: 

Y JJt - (2J'+1) I (2L+1)(J J J') 2 a(L,J) 2 Y. q , J' < J. (2.3) 

The correction factor Is given by 

A(L, J) - [6+a(L) 2 ]/[6+fl(J) 2 ), (2.4) 

with 

fl(n) - 0.13 n B l / (uT) 1 / 2 , (2.5) 

where the rotation constant B Is In cm -1 , the colllslonal reduced mass u is in 
atomic mass units, the kinetic temperature T Is In kelvlns, and the scattering 
length l is In A, with a typical value 1 “ 3A. 

The ylo needed In Equation (2.3) have been obtained from new scattering 
calculations. Using the accurate lntermolecular forces of Thomas et al. 

(1980), we computed cross-sections at eleven collision energies ranging from 
100 to 5000 cm -1 using the IOS scattering approximation, which has been shown 
to be accurate for this system (see Goldflam et al. 1977 and Green 1979). Rate 
constants ylo for L < 32 were obtained by averaging these cross-sections over 
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Boltzmann distributions of collision energies at the appropriate temperatures. 
The resulting values, which are given In Table 1, were then used with Eq. 

(2.3) to provide all the necessary colllslonal rates. (Note that Ylo with L > 
32 were set to zero.) 

It should be noted, finally, that all of these calculations have 
specifically treated excitation of CO by collisions with He atoms. It has 
been argued (e.g. Green and Thaddeus 1976) that excitation by H 2 Is quite 
similar, the main difference being the greater velocity of H 2 at a given 
temperature due to Its smaller mass. This factor will increase the rates for 
H 2 excitation by about 50X from the values for He. Experimental evidence in 
support of this procedure can be found in available pressurs broadening data 
(Nerf and Sonnenberg 1975; see also Brechignac et al 1980). 

The population of each rotational state can now be calculated from the 
calllsional rates Yjj* by solving the equations of of statistical equilibrium 
for these states. For the particular case of shocked CO, certain simplifica- 
tions can be made: 

1) Each line Is assumed to be optically thin. 

2) Infrared pumping by dust emission Is assumed to be negligible. 

3) Only the ground vibrational state contains a significant population. 
Observational evidence to date supports these assumptions. The equations of 
statistical equilibrium are set up as follows. In steady-state, the rate at 
vhich molecules leave each state Is equal to the rate at which they enter 
It: 

n j( A j + n ^ H 2^ j» y JJ») “ ( n ^ H 2^ ^J* J n J* + n J+l A J+1 ) " 0 (2.6) 

where n Is the density of molecules In state J, n(H 2 > Is the density of H 2 
J 

molecules, and Aj Is the A-coefflcient for a radiative transition from state 
J to J-l. There is one such equation for each rotational state J. Together 
with the normalizing condition, J n - n(C0), these expressions make up a 
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complete eet of linearly Independent equations, the solution of which Is 
straightforward. 

In the optically thin case the line emission coefficient Ij for the 
transition J to J - 1 is given by: 

Ij ■ A j nj/n(CO) {erg sec 1 molecule 1 ar *1 (2.7) 

wtv«re v, the transition frequency, Is just 2J times the rotational constant B 
(neglecting centrifugal distortion). Figure 1 shows the line emission coeffi- 
cients as a function of J for J < 50, at several temperatures and densities. 
These curves result from computer solution of the equations (2.6) with J < 60, 
as described above. 


III. ANALYTIC ESTIMATE OF POPULATIONS 

In order to make progress In obtaining an analytic estimate of the level 

populations of the CO molecule, It is necessary to make a drastic 

simplification in the collision cross-sections ojj* connecting levels 

J and J' . Here we follow HM and adopt 

Yjjt - <ojj. v> - gj» exp(~ Eji/kT) ovt ( 3 . 1 ) 

Z 

where Z Is the partition function, vj ■ { 8 kT/mn(H 2 ) ] Is the mean speed of 
the Hi molecules, gj* - 2J' 4 - 1 is the degeneracy of the state J' , and a is 
the total collision cross-section out of level J: 

1 ^°JJ’ V> ■ ovj, (3.2) 

J-0 

Equation (3.1) for the cross-section has the advantages that it satisfies 
detailed balance and it facilitates solution of the level population equations, 
tlons. Its major disadvantage Is that it depends only on the final state 



and therefore omits a conspicuous feature of the calculated cross-sections, 
namely that collisions with small AJ are strongly favored over those with 
large AJ (cf. Table 1). Numerical evaluation of o shows that It la about 
10-15 co 2' ag ai8um ed by HM; the effective value required to determine the 
level population Is smaller, as we shall see below* With expression (3*1) for 
the cross-section, equation (2.3) for the level populations readily simplifies 
to 

nj(Aj + 1 n a vr) - nj+i A J+1 + 1 n nj* o vj (3.3) 

2 2 


where 

nj* - n(C0) SJ e~ E j/ kT (3.4) 

Z 


Is the LTE population, n Is the number density of hydrogen nuclei, and n(H2) * 

n/2 in a fully molecular gas. In deriving equation (3.3), we set Ojj - 0 

although it Is non-zero according to equation (3.1); we also set 

I nj. <• n( CO) (3.5) 

J’*J 


Although in this paper we are focussing on the case of optically thin 
emission, It Is a simple matter to include the effects of line opacity via the 
escape probability formalism (e.g., deJong 1973, HM). If ej < 1 Is the 
probability that a photon In the J ♦ J-l transition can escape the emission 
region, then the transition probability Aj Is replaced by Aj cj In equation 
(3.3). Making this substitution, we Introduce 


aj i 2A J C J (3.6) 

novj 

Then, in terms of the departure coefficients bj - nj/nj*, the population 
equation becomes 


9 


This equation can ba readily solved by changing the dependent variable 
to ajbj. For J > 1, the result is 

- .* J ' 

ajbj ■ l «j' *** *J" (3.8) 

J'-J nj w J"»J 1 + ej" 

The ground state population is then given by 

* 

b o " b l M ill* 1 (3.9) 

V 

In view of the approximate form of the cross-section used in obtaining 
the solution (3.8), considerable simplification is warranted. For simplicity, 
assume J » 1 so that 


where, for CO, Kq 


Aj - Ao j3 
Ej * E e j2 

1.118 x 10" 7 s” 1 and E 0 /k - 2.76b K, 

y = Ej . E 0 J 2 = j2 

kT TT 


(3.10) 

Define y and Jj by: 

(3.11) 


For kT » E Q , the partition function is simply Z - Jj 2 . The energy difference 

AEj between levels J ind J-l is about 2E 0 J so that 

AEr/kT « 2y/J (3.12) 

The quantity aj defined in equation (3.6) can be expressed 

«j - j y ej (3.13) 

4 n 


where 


. 2 . 

n cr = 8 JT Aq (3.14) 

OVj 

is the density of hydrogen nuclei at which colllslonal de-excitation becomes 
Important in an optically thin gas (HM). 

First, consider the high density limit in which aj+0 and the levels are 
nearly in LTE. Treating both aj and J" ^ as small and keeping terms to 2 nd 
order in small quantities, we find that the solution (3.8) reduces to 


1 
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b J ■ 1 + (n cr Cj/n) y 
1 + l(nj r tj/n)y* 
7 


(3.15) 


Next consider the low density Holt In which ij ♦ •», so that 

«jbj* l 

J’^J nj* 


(3.16) 


Replacing the sua by an Integral, we find for J » 1 


b j ^ 

~T (n cr ej/n)y 2 
2 


(3.17) 


A generalization which Is valid it J * 0 as wall Is 

bj ■ 1 

V r + i (ncr e j/n) y 2 
2 


(3.18) 


There Is considerable latitude In combining equations (3.15) and (3.18) 

Into a general approximation for all densities. Since the quantum mechanical 

calculation of the crosa-aectlons shows that o decreases slowly with J, we 

expect that n cr will Increase slowly with y. Hence we make the replacement 

"cr e J ♦ wy* (3.19) 

n 

and treat w and S as parameters to be determined by comparison with the 
numerical results for the level populations presented In S II . A simple 
expression which reduces to equations (3.15) and (3.18) in the appropriate 
limits Is 


b J 


Sr 


>o 1 + 1 wy 

7 




(3.20) 


At high densities, b Q 1, as Indicated In equation (3.15). 

Vith this expression as a foundation, we have developed the following 
approximation for the departure coefficients in the optically thin case by 
trlal-and-error comparison with the detailed numerical results: 
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f » , 1*3 . 

bj - f i + »y_ ] i 

bo“l -f 1. wy2»3 (l + yyl«^) f 
b 0 - (1 +2 w 0.83)l/2 


(3.21) 

(3.22) 


f - (1 + 0«55 y 2 - 2 )2 

T3-0.2 + 6w“°.75 


(3.23) 


w-2 I 75jr^ (3.24) 

n 6 

a - 0.74 + 0.16 T 3 I /2 (3.25) 

where ng = n/(10^ cm"^) an< j = t/( 10^ K). The factor f he* been Inserted to 
reduce the populations at large value* of J, thereby accounting for the In- 
efficiency of collision* with large dJ. There Is an upper bound on b 0 deter- 
mined by the condition that all CO molecules be In the J-0 level, which gives 
b 0 < Z ■ Jt 2 ; however, this limit is reached only for densities below that at 
which the above approximation Is valid (see below). At T 3 ■ 1 and J ® ■ 

19, the value of w In equation (3.24) corresponds to n cr > 2.8 x 10 * cm~3, or 
0 « 4 x 10~16 cm^. As mentioned above, this effective value of the cross-sec- 
tion Is several times smaller than the total de-excltatlon cross-section. 

This analytic approximation is remarkably accurate over a wide range of 
physical conditions. Let Jp be the rotational quantise msaber of the highest 
level which produces a line with an Intensity within a factor 50 of the 
strongest line; If the strongest line Is at J B , then typically J p - 2J B . Com- 
parison with the numerical results shows that our approximation Is correct to 
within a factor 1.33 for 0 < J < J p < 50 over the range 2 x 10^ cm~3 < n and 
250 K < T < 2000 K, except at low densities and high temperatures: at n ■ 

2 x 10 5 cm** 2 , the accuracy drops for T > 1750 K and Is 1.4 at 2000 K; at n ■ 

2 x 10 4 K, the accuracy drops for T > 1000 K and Is 2.0 at 2000 K. 
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The accuracy of th« approximation was checked at 5 densities, n ■ 2 x 10^ 
cm“3 f 2 x 10^ cm"3, •••2 x 10® cm"^, and at lntarvala of 250K from 250 K to 
2000 K. .-\r T - 100 K, tha accuracy la 1.5 for n > 2 x ’0 4 cm" 3 . 

Tha approximation davalopad hara la raadlly axtandad to tha casa of fl- 
nlta optical dapth in tha llnaa by multiplying tha valua of w glvan in aqua- 
tion (3.24) by the aacapa probability cj (aae HM). Ue have not attempted to 
aaaaaa tha accuracy of thia procedure, however. 

IV. POPULATION INVERSION 

Population lnveralona result in maser action and can produce extremely 
intense sources of radiation. Because of high abundance of CO in molecular 
clouds, Inversion of CO levels could be particularly Important. Goldsmith 
(1972) suggested that the J • 1 level of CO could become Inverted because 
collisions from J • 0 to J ■ 2 occur at a rate comparable to those from J ■ 0 
to J “ 1, whereas radiative transitions from J • 2 to J ■ 1 occur much more 
rapidly than from J ■ 1 to J • 0. In a detailed calculation Including 
radiative transfer, Leung and Liszt (1976) found that at 40 K the 1*0 
transition became euprathermal (excitation temperature > kinetic temperature) 
but not Inverted. Because of the lower optical depth in l^CO, much higher 
excitation temperatures were achieved in l^CO than In ^C0. Actual Inversions 
were predicted In the CS molecule at this temperature (Liszt and Leung 1977). 
Calculations of the populatlone In cyanoacetylene (HC 3 N), which Is a linear 
molecule like CO, Indicate that the 1*0 transition can be inverted for 
lO^ cm" ^ < n < 10 ^ cm"3, whereas the 2*1 transition Is inverted over a 
narrower range (Morris et al. 1976). The latter authors claimed evidence for 
a weak cyanoacetylene maser In Sgr B2. 
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Using ths detailed msnerlcal calculations discussed In III, we have nade 
a aysteaatlc search for population Inversions In CO over a wide range of 
temperatures and densities. We have not Included any radiative transfer 
effects, because that would entail making additional assumptions about the 
geometry of the source and the Intensity of the background radiation. We find 
that the 1+0 transition can be Inverted for T ** 50 K over a narrow range of 
densities; as the temperature Is Increased, the range of densities over which 
an Inversion occurs Increases and higher levels become Inverted as well. The 
mechanism for producing the Inversions is a direct generalisation of that 
suggested by Coldsmith (1972): colllslonal rates vary more slowly with J than 
do the radiative rates, so It Is possible for molecules to become trapped at 


high J. 


Our results are stamarlzed In Table 2 and Figure 2. The magnitude of the 
inversion njgj'V(nj.i gj-l~*) and the ‘‘saturation column density" N ga t are 
presented as a function of J, n(H2), and T. When the CO column density equals 
N t at> the optical depth in the J + J-l transition is *1. For N(C0) < N gat , 
our optically thin approximation Is valid but the amplification Is small. For 
N(CO) > N ggt , the amplification could be large In principle, but a detailed 
radiative transfer calculation would be needed to make certain. The value of 
N ggt is determined by the optical depth tj In the J ♦ J-l transition. At line 
center we have 


T 


J 


N 


J-l 


>J * 3 *J 
gj-1 8w 3 f £ tv 


(1 - N J «J-I ) 

Nj-lj 


(4.1) 


where Nj Is the column density In level J. The line width Av is taKen to be 
the half width to 1/e points of the line, expressed In velocity units; for a 
thermal profile, Av - 12.9 (T/A) 1 ^ 2 la s**l, where A Is the molecular weight. 
Noting that Aj « y 2 X* 3 , where y la the dipole moment, we can solve equation 
(4.1) for N(C0) when t • -1 and obtain 
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N „ - 2.74 x 10 13 [ N J-1 J ( N J 8J ~ 1 - l)] -1 cm* 2 (4.2) 

n7C0T (2J-1)/T V Nj-! gj 

where the ratios Nj_j/N(CO) and Nj/Nj_i may be obtained from the calculations 
described in $11. The values for N aat given in Table 2 are minimum values 
because we adopted the minimum line width - the thermal width. CO has 
unusually large values of N gat because of its small dipole moment. 

The conditions required for a population Inversion may be understood as 
fellows. Let J £ be the level at which nj is a maximum; obviously, an 
Inversion is possible only for J < J*. If the density is too high, then the 
populations approach their LTE values and no inversion is possible. If either 
the density or the temperature is too low, then J* will drop below J and again 
no inversion is possible for level J. The upper limit on the density can be 
estimated from the approximation developed in Sill; unfortunately, the lower 
limit occurs at a density too low for the approximation to be valid. Since an 
inversion is possible only for low values of J, the quantity y ■ J 2 /J>r 2 is 
small and the factor f can be set equal to unity in equation (3.21). Hence bj 
depends only on the two variables w and y, and not on n, T, and J separately. 
The condition that an inversion occur for a specified y can than depend only 
on w. Since an inversion is possible only for densities below some critical 
value (which depends on T), the condition must be a lower limit on w; l.e., we 
require w > w cr ^ t for an Inversion. Inspection of Table 2 suggests w cr ^ t » 4, 
which gives 

n(H 2 ) < 10 5 * 5 T 3 a cm -3 (4.3) 

as the condition for an inversion. At T ~ 50 K, the upper and lower limits on 
the density converge at n(H 2 ) ~ 3 x 10 4 cm -3 ; the numerical calculations show 
that an Inversion is impossible for T < 50 K. As T increases above 50 K, J* 
rises so that higher levels become inverted as well. 


IS 


As emphasized above, our calculation of population Inversions Is valid 
only in the optically thin limit, N < N gat , which seemingly precludes large 
amplifications. However, if the CO Is in a thin sheet, then our calculations 
may be directly applicable even if the amplification is large: in most 

directions one could have N(CO) < N sat so that the level populations are 
unaffected by radiative transfer, but when viewed edge-on the sheet could have 
N(CO) » N sat , producing a source of high brightness. The optically thin limit 
may also be preserved if the gas has a large superthermal internal velocity 
dispersion. Such conditions are expected for interstellar shocks. 

Observations of regions containing Interstellar shocks, with sufficient spatial 
resolution to resolve the shock thickness, could reveal regions of anomalously 
high brightness in the lower CO rotational lines. A possible example of such a 
region is IC 443 (De Noyer 1979a, b) where some components may be edge-on, 
although the excitation seems too small to produce inversion. Unfortunately, 
in Orion, the best-studied example of interstellar molecular shocks, the 
observed CO emission originates in gas which is too dense (n[H2l - 10 6 cm -3 ) to 
produce a population inversion. 

Detailed calculations are required to determine if significant amplifica- 
tion can occur in actual Interstellar shocks. Reference to Table 2 shows that 
relatively large column densities of warm gas are required to produce one e- 
foldlng in the intensity: N(H) - IQ 20-21 cm" 2 at n(H 2 ) ■ 10 4 cm" 3 and N(H) - 
10 21,5 " 22,5 cm" 2 at n(H2) * 10 3 cm" 3 , under the assumption that n(C0)/n ~ 10" 4 . 
If the emission occurs in a magnetic precursor to the shock (Draine 1980) then 
n is the pre-shock density; if the emission occurs in the shocked gas, then n 
1 s the shocked density, which is typically 10-30 times greater (HM). 
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V. OPTICALLY THIN COOLING RATES 

Because of Its relatively high abundance, CO Is often an Important cool- 
ant In molecular clouds (e.g., Dalgarno, et al 1975) • We define the rota- 
tional cooling rate coefficient L roC by setting n n(C0) L rot equal to the to- 
tal CO rotational emission rate; assuming that the gas is optically thin, we 
have 

n n(CO) L rot « £ nj Aj AEj - 4im(C0) l Ij (4.1) 

J-l J-l 

where AEj is the energy difference between levels J and J-l and Ij is defined 
in equation (2.7). HM obtained the approximate result 

L ro t • 4(kT) Aq (4.2a) 

"Ml + (°cr/ n ) + 1.5( n cr/ n )l/ 2 ] 

_1 

“ 2 k T ^T (4.2b) 

1 ♦ ("/ ncr) + 1.5 (n/n cr )l /2 

where the equivalence of the two expressions follows from the definition of 
n cr , equation (3.14). The factor 1/2 in equation (4.2b) is simply n(H 2 )/n. 

Using the numerical results discussed in III, we find that this 
expression for the CO cooling rate is quite accurate. Setting 

n cr - 3.5 x 10 6 T 3 2 / 3 cm” 3 , (4.3) 

which corresponds to o - 1.6 x 10 " 16 cm 2 at T 3 - 1, we find an accuracy of 20% 
for 2 x 10 * cm " 3 < n < 2 x 10® cm " 3 and 250 K < T < 2000 K. However, HM's 
numerical value for the critical density (n cr * 6 x 10 3 T 3 I / 2 cm" 3 ) is almost 
an order of magnitude too low, primarily because they set 0 equal to the total 
de-excitation cross-section of about 10"^ 3 cm 2 ; in addition, thi*y used a ther- 
mal velocity appropriate for atomic hydrogen. This underestimate of n cr would 
lead to an overestimate of the CO cooling at low densities by almost an order 


of magnitude 
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In order to see how the o In equations (3.14) and (4.2b) Is related to 
the colllslonal rates presented In Table 1, we proceed as follows. In the 
low-density limit 


^ 1 he *. •' Ej/kT *, 

n J J® J J 

(4.4) 

/ y J0 yl* 

“ “A 

(4.5) 


For 2 < J < Jf, the results In Table 1 can be approximated as yjo 8 Y20(2/J) 2 ; 
for J ■ 1 and J > J*p, Yjo Is smaller. Evaluation of the Integral gives L r0 £ - 
2 kTT 2 ()C, where the factor C allows for the smaller values of ijq for J > 

(y > 1); we somewhat arbitrarily set C « (1 - e“ l ) » 0.6. Comparing this 
result with equation (4.2b) in the limit n << n cr , we obtain 

ov-p a 2.4 120 (4.6) 

If allowance is made for the mass difference between He and K 2 , one finds that 
120 " 3.2 x I0 " 11 T 3 1 / 3 cm^ 8 -l f or gj. Finally. Inserting these results Into 
the deflrltion (3.1 ) of n cr gives n^ ■ 4.3 x 10 6 T 3 2 / 3 cm -3 , quite close to 
the value found above. Equation (4.6) should allow one to estimate n cr for 
molecules other than CO. provided their excitation rates behave similarly. 

Finally, we note that equation (4.5) applies only In a gas of molecular 
hydrogen. Behind fast shocks, relatively large concentrations of atomic 
hydrogen can build up (Hollenbach and McKee 1980). If X 2 • n(H 2 )/n, then 
In the presence of atomic hydrogen the critical density (eq. 3.14) becomes 


l cr 


4 J 2 A 
T 0 


^T tx2 + ( 1 “ 2x 2) 


(4.7) 


there or, Is the effective cross-section for rotational excitation in 
collisions of atomic hydrogen with CO. In a fully molecular gas, X 2 ■ 1 and 

7 

equation (3.14) Is recovered. 
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VI. CONCLUSION 

The msnerical and analytic result* we have presented for the rotational 
spectrum of CO should be useful in interpreting observations of Interstellar 
CO emission from warm molecular gas and in generating theoretical models of 
shocks in molecular clouds. The results are valid for optically thin CO which 
is colllslonally excited by H 2 collisions; infrared pumping has been ignored. 
The analytic results can be extended to the optically thick case in an approx- 
imate fashion by using the escape probability formalism. 

An important result of our analysis is that population Inversions in the 
lower rotational levels of CO can be expected to occur in warm molecular gas 
with density in the range n(H 2 ) ~ 10 3 ” 3 cm - ^. The number of levels for which 
an inversion occurs the column density required for significant amplification 
both Increase with temperature (see Table 2 and Figure 8). Interstellar 
shocks appear to be promising as sites of CO masers. The expression for the 
CO cooling rate obtained by HM has been verified. However, we find that the 
effective cross-section which determines the cooling rate is significantly 
smaller than they assumed (o - 4 x lO - "^ cm^ rather than 10"*5 CB 2) f and that 
the cooling rate at low densities is correspondingly reduced. The effective 
excitation cross-sections for OH and H 2 O, the other dipolar molecules which 
are Important coolants, are probably also less than the value 10~15 cm 2 
assumed by HM; however, since the AJ required to obtain a given energy differ- 
ence is much smaller for these molecules than for CO, we anticipate their 
effective cross-sections are closer to 10~15 cm 2 chan that of CO. 
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coefficients for He-CO collisions are given. H2-CO rate coefficients are obtained by Multiplying 
:S in this table by 1.37. 
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Figure Captions 



Lins mission coefflcltnt Ij for ths CO rotational transitions 
J * J-l, as a function of J for sever*! tanparaturas and molecular 


hydrogen dansltlas n(H2>« Tha Ij distributions cobs fron computer 
■ol'lon of tha statistical equilibrium aquations, as discussed in 


III. Also shown, for comparison, ara tha Ij distributions resulting 
from thermal equilibrium (LTE) level populations. 


Figure 2 Population inversion njgj.^/nj-igj as a function of temperature and 
molecular hydrogen density for the J ■ 1 and J ■ 2 states of CO. 
These curves were generated using computer solutions of the 
statistical equilibrium equations, as discussed in III. 
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